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Abstract 

I propose an approximation scheme for asymptoticahy free field 
theories combining both weak couphng and strong couphng series. 
The weak couphng expansion is used to integrate the high frequency 
modes and the resuhing low energy effective theory is solved using the 
strong coupling expansion. In some models there exists an intermedi- 
ate scale at which both expansions make sense. The method is tested 
on few low dimensional models for which an exact solution is known. 
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In strongly coupled theories like QCD there are plenty of qualitative phe- 
nomenological models (quark, bag, string models, strong couphng expansion 
...), but just one approximation which is under control: the asymptotic weak 
coupling perturbation theory. This is only applicable at energies much higher 
than Aqcd- It is commonly believed that for low energy physical quantities, 
like hadron masses, there is no small parameter in which one can expand. 
The theory apparently does not contain a natural small parameter and the 
expansion in a reasonable small "coupling" fails at small energies. 

Nevertheless, the relative success of one of the "phenomenological" ap- 
proaches seems to indicate that it contains at least part of the truth. This is 
the strong coupling expansion which is usually performed using lattice reg- 
ularization with lattice spacing ~ 1/ (few Aqcd) Confinement and chiral 
symmetry breaking are easily seen and the intricate spectrum is nicely ac- 
counted for, at least as far as ordering and quantum numbers are concerned, 
and even the mass ratios are reasonable 0. At one time there was a hope 
that one can extrapolate from the strong coupling series into the "scaling 
region" using Pade approximants 0. This would convert the method from 
just a phenomenological model to a systematically improvable scheme. Un- 
fortunately, this method includes too much guesswork to be controllable and 
it was impossible to reach the scaling region from the apparent radius of 
convergence of these series (about /? ~ 5 in SU{3) Yang - Mills 0). As 
it is usually used, the strong coupling series have very little to do with the 
(asymptotically free) continuum limit we are interested in. Therefore one 
should consider it just as a model for low energy phenomena with the UV 
cutoff M and the coupling g which should then be related to Aqcd- 

Even in much simpler asymptotically free theories the situation is similar. 
In the d=2 a - model and the Gross - Neveu model strong coupling expan- 
sions were performed |^ and qualitatively reproduce quite nicely the exactly 
known low energy features. Yet the continuum limit is still unreachable. 
The available exact S - matrices in these models show what kind of non- 
analyticity in coupling constant is expected. The general dependence of a 
physical quantity, for example a S - matrix element, on a physical coupling 
g has the following structure 0: 

oo oo 

fig) = E E «.n,n^?™e-"/^ (1) 

m=0 n=0 

where a^.n are numerical coefficients. It contains both polynomial and non- 
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analytic in g (but analytic in 1/g) "instantonic" exponential factors. If one 
would be able to separate these dependencies, by expressing / as 



so that -F(a, (3) has a viable expansion in both variables, this would provide a 
systematic method to approximate these theories. Here the "weak coupling" 
a{g) is analytic in g and the "strong coupling" I3{g) is analytic in 1/g. 

In this paper I formulate such an attempt and test it on low dimensional 
QFT like the just mentioned two. Of course, one can wonder whether this 
is at all possible. Apriori it looks like the domain of applicability of the 
expansion in (3 has a very little chance to intersect with that of the expansion 
in a. Indeed, if g is small, 1/g is large and vice versa. Fortunately, this 
does not always imply that for small a{g), j3{g) should be large and vice 
versa. Just a glance at strong coupling expansion results [0] shows that, 
although the asymptotic scaling region is out of reach, the effective weak 
coupling near the strong coupling radius of convergence is quite small. In the 
SU(?>) YM theory the radius of convergence [3max = ^/Qmin ~ ^ corresponds 
(using naive perturbative RG) to a relatively small effective weak coupling 
expansion parameter a = g'^/4:TT ~ 0.1. Even taking into account the fact 
that the a is a lattice one (not the MS), this corresponds to the values 
at which perturbation theory for high momentum quantities is supposed to 
work. The same happens in two dimensional asymptotically free theories. 
Consequently, there is a (albeit small) window in the coupling at which both 
a and (3 are small enough to produce a reasonable series. The "loop factors" 
l/(47r)^ in the practical weak coupling expansion parameter a{g) are partly 
responsible for this. But this should be checked in any particular case since, 
for example, the symmetry group factors tend to reduce the windowQ. 

Now, I will outline the general method to calculate a physical quantity. 
The model's lagrangian on the UV cutoff A scale is proportional to a local 
dimension d operator: C\{(f)) = -^Co{(p), where g^ is a bare coupling con- 
stant and (f) are fields^. 



^The location of the "window" and even its existence may depend on the choice of 
coupHng and therefore on the regularization and the normahzation scheme. 

^ We consider here only models like pure Yang - Mills theory or chirally invariant QCD 
without additional couplings. The generalization to more complicated asymptotically free 
theories is trivial. 
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1. Aa a first step, one integrates the high frequency modes from A down 
to a certain scale M. M should be larger than Kmodei , so that the weak 
coupling expansion (or the loop expansion) is applicable where Amodei is the 
corresponding dynamically generated scale of the "model". The result is, of 
course, a complicated functional: an effective action of the type extensively 
studied recently 0. It typically allows the standard derivative expansion in 
powers of of the following form: 

oo oo -1 

A/[0]= E E%^?\t7^A(0) (3) 

i=-lj=0 

where Cj is a dimension d + 2j local operator. The first operator Cq is, of 
course, the original lagrangian now containing however the low frequencies 
only. If the regularization is the lattice one, it is simply a coarse grained 
version of the original lagrangian. Its coefficient is a scale dependent coupling 
at the scale M: ^ = ^+I]i=o ^io9\ can be used as an effective expansion 
parameter. One then can reexpress g\ via qm in the other terms of eq.(|^) as 
well and take the continuum limit (7a ^ 0, A ^ oo with Amodei fixed. Now 
we introduce j3 = I/qm and a which is analytic function of qm (typically 
a small coefficient times qm or g^). Then the effective lagrangian takes a 
form: 



Cm\(P]=P. 



M 



oo oo 

^0 + E E (^ijOtMTHJ^ji(P) 

i=l j=l 



(4) 



An arbitrary physical quantity, like an expectation value of an operator 
0(0), can be expanded in (practically, in corresponding um) and in 
(which means in (relevant momentum)^/M^) as: 

r 0((h)e~-f oo oo 1 

<0>= =< O >, +/3mEE%«mT7^ < oc,{<p) >, +... 

(5) 

where < ... >c means the connected part of the expectation value in the 
coarse grained version of the original theory and the higher order terms come 
from the expansion of the exponent. 

2. Now the strong coupling expansion is applied to calculate these VEVs 
assuming that the strong coupling expansion parameter j3 = I/qm is inde- 
pendent of a in eq.(|D. Note, that the leading order term in a coincides with 
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the conventional " phenomenological" strong coupling model, in which the 
inverse lattice spacing M, however, is limited to values inside the strong - 
weak applicability window. 

The complexity of such a calculation depends on the quantity and the 
precision one would like to achieve. Usually strong coupling series are easy 
to evaluate to a very high orders. Consequently the scale M can be chosen 
in such a way that (3m is just below the strong coupling expansion radius 
of convergence. Simultaneously, the scale M should be sufficiently high or 
alternatively the relevant energy scale sufficiently low so that just a few orders 
in derivative expansion are needed to achieve the desired accuracy. This 
means that the method is limited to low energy quantities only. The loop 
expansion for the effective action coefficients is probably the most tedious 
part. An example of such a complexity estimate is given later. 

Rather then describe the general method in more detail, I will proceed to 
apply it to a few simple models in which it is very transparent. I will start 
with the simplest possible one: the Ising chain 



considered as a ci = 1 QFT. Asymptotic freedom in language of statistical 
mechanics means phase transition at zero temperature. This occurs in d=l 
Ising model and generally in d=l 0{N) symmetric fixed length spin systems. 
In d=l nearest neighbours interaction spin systems the RG transformation 
can be performed exactly and then compared with the weak coupling (low 
temperature) expansion. This generally is not available even for d=2 models 
for which the S - matrix is known. 

The exact " decimation" of the Ising chain from the scale a = 1/A to A = 
1/M is especially simple. Not only does the interaction remain the nearest 
neighbours one, in addition the coarse grained effective action contains just 
the original term with a new coefficient H: 



— = Arctanh 
9m 



Tanh^^/^ ( - 



1 - 1 Log (^] + l^e-/- + O (^e-^/A (7) 
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The perturbative expression is obtained by neglecting exponentially small 
pieces: I/qm = 1/5'a — iLog(A^/M^). This is the usual logarithmic scale 
dependence, similar to that in the d=4 YM theories or the d=2 Gross - Neveu 
and a models, corresponding to the following /3 function: P{g) = —\g'^- Note 
that there are no contributions to the (5 function higher then the second 
order in the coupling. The corresponding A parameter A/ is defined by: 

A/ = Ae ''(f) = Ae '^'^ The perturbative expression takes into account 
only the lowest energy configuration, Fig.4a,b, while, as will be shown later, 
the exponentially small pieces account for kink - antikink pairs. Fig 4c, d. 
Since the effective action contains only one operator rather than the usual 
infinity of various operators, no derivative expansion is needed. Now to 
complete the step 1, we just have to express qm via (^a and A/ and take the 
limit A ^ cx). The result is: I/qm = |Log (M'^/Aj). 

To perform step 2 for the mass of the single existing excitation one sim- 
ply exactly solves the effective model and then uses the exact expression to 
generate the strong coupling series in Pm = ^/qm- 

m/M = -Log[Tanh(/3M)] = -LogiM + \f3lj ' I^m + ^Pli + - (8) 
The deviation of this from the exact value, mexact/Aj = 2, for various M is: 



A, 3VM/ nlij ^'^lUjj <''^ 

The radius of convergence of the strong coupling series in eq.(§) is Pmax = 
tt/2. At this point Mmax = e'^Aj ~ 23. lA/. The maximal precision one can 
achieve is therefore {niexact — ^)/^exact ~ 10^^. 

In this very simple example only strong coupling expansion came into 
full play. The other two ingredients, namely the weak coupling expansion 
and the derivative expansion, have not. The weak coupling expansion for 
I/qm terminated at the one loop level and higher loop effects for physical 
quantities are just the RG improvement, while the derivative expansion was 
not needed. 

Now we consider a slightly more complicated solvable model in which 
these two elements already appear: the 0(3) symmetric spin chain: 



Za= ei^^-^^^^+^ (10) 
52=1 
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This model has a phase transition at zero temperature, although the asymp- 
totic freedom is powerwise rather then logarithmic: f3{g) = —g + ... [^]. 
Using the transfer matrix method, one is still able to perform an exact dec- 

imation. The infinite dimensional transfer matrix K{S, S') = can 
be diagonalized in the basis of spherical harmonics |H|: 



1=0 



Im \ 



T^gA 



1=0 1=0 

Fi are eigenvalues of the matrix. Here 6 and are the spherical coordinates 
and /z+i/2 and Pi are modified Bessel functions of the first kind and the 
Legendre polinomials respectively. 

After the decimation the effective action still contains only nearest neigh- 
bours interactions, but now all the powers of ^DS*, where □ is the lattice 
second derivative, are present: Cm = const + j^SOS + J2i=i Ci{SOSy. The 
coarse grained lattice's transfer matrix is given by the (A/M)s power of K\. 
It means that the eigenvalues are: F/^ = {Fj^)^/^ . Introducing Kh = A/gA 
and taking the continuum limit, one finds the effective action 

F,^ = e-^^. (12) 
As in the case of Ising , this in fact defines a perfect action |TT[. Generally, 



for a nearest neighbour model with eigenvalues Fi the spectrum of excitations 
is given by 

mi/M = Log . (13) 

Using this one can easily check that all the continuum masses, which are 
simply levels of the quantum rotator, are correct: rrii = ^^^^Ajj. The coeffi- 
cients of the first few terms in the derivative expansion of an action defined 
by the transfer matrix F are determined by following equations: 

oo 

z = E(2/ + l)Fi 

1=0 

z/gM = f:i2l + l)^-^^Fi 

1=0 ^ 
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where z is a physically irrelevant constant. Taking here the action defined 
by eq.(|12|) and using asymptotic form of the Bessel functions [|I0|, one finds 
that the original SOS dimension one term comes with a coefficient I/qm = 
1/^-1/6-^/180, while that of the next one is: ci = -1/3^+7/90+11^/2700. 
Here g = Ah/M. 

Now the resulting model can be solved exactly. To leading order in the 
derivativeand the weak coupling expansions, one gets the known result for the 
energy gap 0: m = — MLog [Coth(l/^) — g]. This can be used to generate 
the strong coupling series. This should be compared with the exact value 
rriexact/ Ah = 1- The difference expanded in inverse powers of M is: 



A 



H 




(15) 



The strong coupling series radius of convergence is now a bit larger than 
that for the Ising chain, Pmax = ^r. The corresponding maximal scale is 
however very low, M^ax = ttA//, due to the fast running of the coupling. In 
order to obtain better precision, more orders in the expansion are needed. A 
peculiarity of this model is that weak coupling corrections are of the same 
order as the derivative expansion corrections since the asymptotic freedom 
is a powerwise one. So they should be computed together. I calculate here 
only the next to leading correction. 

For this, one should take into account the next to leading order con- 
tribution in g to l/gM and the leading contribution to the coefficient of 
the four derivative term. To find the spectrum one should reconstruct the 
eigenvalues Fi from which the mass is obtained using eg. ([131). This can 
be done using the completeness relation of the Legendre polynomials [|I^]: 
Fi = J^^dye^^^'^^^y~^^'^'^'^^y~^^ Piiv)- To current order this integral still can be 
expressed in terms of the error function, but it is more instructive to perform 
it perturbatively in g: 



^0 = J ^ c?l/e 



^1 = y ^ dyye 



bg 

i-ky-i)-^{y-ir 

D bg 



b 



g + ]:~g' + ... (16) 

6 
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Using the formula for the energy gap eq.(|13]), one gets the correct result 
within a precision of ~ 0.1 and the first term in eq.(^) disappears. 

This can be continued to higher orders. 

Now I turn to a two dimensional logarithmically asymptotically free mod- 
els which have more similarity to four dimensional gauge theories: the Gross 
- Neveu and the a - models. I will describe here the former only, since the 
results are similar (which should not be too surprising, since in many respects 
the models behave similarly). The model defined by the lanrangian 

H = itp^^^'^ + -^{VVf (17) 

(a = 1, ...N) has a rich particle spectrum and exhibits dynamical (discrete) 
chiral symmetry breaking. Although the exact S - matrix of the theory is 
known [0, in order to compare the expansions with exact results at every 
step I will consider the model in the large number of flavours limit only. This 
time, however, the 1/M^ expansion to an arbitrary order will be calculated. 
The first step is the perturbative integration of the high frequency modes. It 
will become clear shortly that only effective potential 

y(^^) = c^M' (^— v^^j (18) 

will be needed to calculate the fermion's mass in the large limit (of course, 
even within the large N approximation there might appear additional irrel- 
evant operator, those containing derivatives, like tpi^Otp or tptpniptp) . Per- 
turbatively, only the one particle irreducible diagrams shown on Fig.l con- 
tribute. 

The internal loops carrying momenta between M and A only are inte- 
grated and external momenta are all zero^ Other IPI diagrams vanish due 
to (unbroken at this stage) chiral symmetry or are of lower order in 1/A^. The 
first, the third and similar chain diagrams determine the relevant dimension 
two term renormalization: 

1 1 1 . f 



9m 




(19) 



■^We use the sharp momentum cutoff regularization which at this level does not lead 
to any ambiguities although the calculation (together with analogous calculation in the 
nonlinear a - model) can be done, for example, on the lattice. Although this is very 
helpful to actually perform strong coupling expansion for finite N p^ , it unnecessarily 
complicates the large N calculation. 
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while the second is the coefficient of the dimension four irrelevant operator: 



<===9J^(5^-^) (20) 



57r 



Reexpressed via Aqn = Ae '"^^'^ and at the continuum limit A — > oo, they 
become: Ci = 5fA//2 and C2 = glf^^. 

At ffist, I neglect the corrections. The mass should be evaluated 

using the strong coupling expansion. Instead, I will find it exactly and then 
expand in Pm- The large N dynamically generated mass mo is the solution 
of the gap equation: 



1 1 pdp 1 f + ml 



Qm tt Jo + rriQ 2n \ tjiq 
In terms of Kgn it is 



(21) 




1 + +••• (22) 



kl^ 1 - {Kgn/MY M2 

and should be compared with the exact value rriexact/ ^gn = 1- 

Now I will show that the leading correction eliminates the sec- 

ond term in eq.p^) thus improving the precision to 0(1/M^). The fermion 
propagator to this order in is 

G{p) = Goip) - < i^ipm-p) jii'i^Y >c (23) 

where < ... >c denotes connected evaluated contributions shown on 
Fig.2. 

All the diagrams (for the Feynman rules, see, for example, [13] 



are 



evaluated with mass rriQ. Note also that the other dimension four operators 
like i/jc^OiIj do not contribute; they all vanish in the large N limit (see Fig.l). 
The result for the mass is (there is no wave function renormalization at large 
N) ml/AQi^ = 1 + 0{1/M^), as advertized. Higher orders in derivative 
expansion gradually help to recover the precise result. 

The radius of convergence of the strong coupling series is evident from the 
exact solution eq.(|2l]) for = ^rw;^: Pmax = l/s'mm = 1- Consequently, 
the maximal possible value of the scale M is then M^ax = ^^-^gn ~ SS.IActv- 
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Now we would like to estimate the practicaly achievable precision of such a 
scheme. Suppose we would like to calculate to order Us in strong coupling 
(so that we neglect the order /J^^^"*"^)), to order in weak coupling and 
to order 2nct in derivative expansion. Then the precision of the mass, A = 
(m - inexact) /friexact, is given by 

A = mmM[A(M)] 



A(M) 



max 



/3 



{ris+l) 



a 



(ra„+l) 



2 \ (rid+l)' 




(24) 



For example, if Ug = 10, = 1 and = 1 one obtains A = 0.045 for 
M = 10.7 Agn and /? = 0.75. If, instead of = 1, we are able to calculate 
riw = 2 (no change in Ug and n^), we improve to achieve a precision of 
A = 0.015 at /3 = 0.68, M = SAAqn- In the first example out of three 
quantities in eq. (p^ the weak and the strong coupling expansion determined 
the optimal M, while derivative expansion was not critical. In the second 
example the derivative and the strong coupling determined the optimal M, 
while weak coupling was not critical. 

The reader might wonder about the following dynamical question. What 
happened to spontaneous chiral symmetry breaking during this calculation? 
On the one hand, we used the weak coupling perturbation theory that does 
not know about the symmetry breaking to calculate the effective potential 
and only later introduced the strong coupling expansion. The strong coupling 
expansion (actually an exact solution of the perturbative effective action) re- 
covered the full effect of the chiral symmetry breaking. We would expect 
that, at least some imprecision would be introduced by the first step. What 
happened here in the large N limit is that the effects of the dynamical sym- 
metry breaking are felt by the effective potential only when the scale M drops 
below Agn- This question is interesting in its own right, so I calculated the 



exact large effective potential using methods first applied in to bosonic 
models. 

Here only the result is presented: 



2{A' 



Agn, 



1 (M^ \ 
-A- -An (25) 
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where u denotes and w{u) is a solution of 



For few values of M it is given on Fig. 3. 

Differentiation the potential at the origin for M > Aqn, one recovers the 
coefficients q to any order. One can see in Fig. 3 how the shape of the 
effective potential changes as more modes are being integrated out. Only 
when the modes below dynamically generated mass are integrated out will 
there be nonperturbative effects. For example, the slope at the origin of the 
potential for all the values M > A^tv is exactly zero, but below Aqn is starts 
rising till it reaches the value of the dynamically generated fermion mass at 
M = 0. 

Let us return now to the double expansion. Notice that, since even at high 
energies there are nonperturbative effects, the use of the perturbation theory, 
even just for integration of high energy modes, should generally introduce a 
systematic error. We have seen this in the Ising chain case. What is the 
nature of these errors and can they be reduced? At least in the case of the 
Ising chain, the answer is that they can be systematically removed if one 
takes into account instantons. Indeed, I will show now that the difference 
between the perturbative effective action and the exact one eq.(|^ is caused 
by the multi instanton effects. The Ising chain is a limiting version of the 
quantum mechanical double well potential. So it comes as no surprise that 
they appear [pTsH . 



Suppose we are decimating from some small scale a = 1/A to a much 
larger scale A = 1/M = ba. When the two adjacent spins on the M lattice 
are parallel, there is just one configuration. Fig. 4a, of minimal energy. 

1 ^ ^ 

Therefore at low temperature, the sum J2si=±i ■■■J2sK=±i 

e^^- " can 

be approximated by Z^^ ~ ca. If the spins are antiparallel. Fig. 4b, the 
lowest energy state with one flip is b times degenerate and the sum within the 

"steepest descent" is Z^^ ~ 6e . The coefficient in the effective lagrangian 
therefore is given by e^/^*^ = Z'^^/Z^^ and coincides with the first two terms 
in eq.(0). 

Now let us take into account one flip - antifiip. For the parallel spins. Fig 
4c, one adds Z^i^ ~ b /2e "a ^ while for the antiparallel spins. Fig. 4d, the 
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(ft-6) 

contribution is Z^^ ~ b^/Qe . Now the third term in the expansion of ef- 
fective couphng in small exponentials of the bare coupling, eq. (|^) is correctly 
reproduced. Therefore if one is to improve on the precision achieved earlier, 
one simply has to take instantons into account. Note, that the instanton 
calculus is performed for small coupling and no problem with overlapping 
instantons should arise. Logarithmically asymptotically free models are spe- 
cial in that contributions of instantons to the effective action is actually of 
the same order in as the perturbative contributions from irrelevant 

operators. Indeed, the exponential factors in effective action in terms of A/ 
are simply inverse powers of M. Consecuently, the precision obtained earlier 
m^/M^^^ ~ 10~^ can be pushed lower by including one kink - antikink effects 
to mVMl, ~ lO-*^. 

In other logarithmically asymptotically free models like the a - model 
or (i = 4 nonabelian gauge theories similar contributions to the derivative 
expansion are expected. They might not show up though in the large N 
limit, as the case of the Gross - Neveu and the a models show. Currently 
this contribution in the 0(3) symmetric a model is being investigated [ll2[ . 

To summarize, the double strong - weak expansion method is proposed for 
a quantitative study of low energy phenomena in asymptotically free theories. 
One integrates the high frequency modes using weak coupling expansion (and 
to further increase the precision, instanton calculus) and then the effective 
theory is treated using the strong coupling expansion. It was applied to some 
solvable low dimensional asymptotically free theories and is able to reproduce 
the exact results. The application of the method to realistic four dimensional 
gauge theories should be addressed next. 

Discussions with T. Bhattacharya, A. Kovner, V. Kushnir, A. Speliotopou- 
los and V. Zhytnikov are greatly appreciated. I acknowledge the support of 
National Science Council of ROC, grant NSC-83-0208-M-001-011. 
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